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ABSTRACT 

An outburst of the accreting X-ray millisecond pulsar SAX J1808.4— 3658 in October- 
November 2002 was followed by the Rossi X-ray Timing Explorer for more than a month. 
A detailed analysis of this unprecedented data set is presented. For the first time, we demon- 
strate how the area covered by the hotspot at the neutron star surface is decreasing in the 
course of the outburst together with the reflection amplitude. These trends are in agreement 
with the natural scenario, where the disc inner edge is receding from the neutron star as the 
mass accretion rate drops. These findings are further supported by the variations of the pulse 
profiles, which clearly show the presence of the secondary maximum at the late stages of the 
outburst after October 29. This fact can be interpreted as the disc receding sufficiently far 
from the neutron star to open the view of the lower magnetic pole. In that case, the disc inner 
radius can be estimated. Assuming that disc is truncated at the Alfven radius, we constrain 
the stellar magnetic moment to /i = (9 ± 5) x 10^^ G cm^^, which corresponds to the surface 
field of about 10® G. On the other hand, using the magnetic moment recently obtained from 
the observed pulsar spin-down rate we show that the disc edge has to be within factor of two 
of the Alfven radius, putting interesting constraints on the models of the disc-magnetosphere 
interaction. We also demonstrate that the sharp changes in the phase of the fundamental are 
intimately related to the variations of the pulse profile, which we associate with the varying 
obscuration of the antipodal spot. Using the phase-resolved spectra, we further argue that the 
strong dependence of the pulse profiles on photon energy and the observed soft time lags result 
from the different phase dependence of the normalizations of the two spectral components, 
the blackbody and the Comptonized tail, being consistent with the model, where these com- 
ponents have significantly different angular emission patterns. The pulse profile amplitude 
allows us to estimate the colatitude of the hotspot centroid to be ^ 4°-10°. 

Key words: accretion, accretion discs - methods: data analysis - pulsars: individual: SAX 
J1808.4— 3658 - stars: neutron - X-rays: binaries 



1 INTRODUCTION 

SAX J1808.4-3658 (hereafter J1808) is the first detected 
accre ti on-powered millisecond pulsar (AMSP) din 't Zand et all 
1 19981 : IWiinands & van der Klij Il998h . J 1808 experiences out- 
bursts lasting a few weeks roughly once in two years, dur- 
ing which the coherent ~401 Hz pulsations are observ- 
able. The neu tron star accretes matter fro m a ~ O.O5i\f0 
brown dwarf I Bildsten & Chakrabartvl 1200 ll) at a 2 hr orbit 
dChakrabartv & MorganI 1998h . T he magnetic field strength is es- 
timated to be about 10^-10^ G jPsaltis & Chakrab artv 1999). A 
recent analysis of the pulsar spin-down between the outbursts 
gives a narrower range of B ~(0.4-1.5)x 10** G CHartman et al, 
l2008al , hereafter H08). Among eight AMSP discovered since 1998 
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showing coher ent pulsations over extended inte rvals during their 
outbursts (se e IWiinandsl |200^ : IPoutanenI |2006| for reviews anc 
iKrimm et alJbOOT *). J1808 is the best studied thanks to its five out- 
bursts and their good coverage with the Rossi X-ray Timing Ex- 
plorer (RXTE). 

The analysis of the broad-band spectra of J 1808 re- 
veals the presence of at least two major components; soft, 
blackbody-like emission belo w 7 keV and a powerlaw tail 
with the cutoff at ~100 ke V jGierlinski. Done & Barrel |2002| ; 
IPoutanen & Gierlinskil |20()3| , hereafter PG03). Both components 
are pulsating at the pulsar frequency and therefore are associated 
with the impact of the accretion stream to the neutron star surface. 
The hard powerlaw is most probably produced in the accretion 
shock, while the blackbody might be the heated neutron star sur- 
face u nderneath and around this shock (see iGierlinski & PoutanenI 
1 20051 for the geometry and observational signatures support- 
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Figure 1. The light curve of SAX J1808.4— 3658 during the 2002 outburst. The flux is computed in the 3-20 IceV energy band. We divide the outburst into 
4 stages: pealc (P), slow decay (SD), rapid drop (RD) and flaring tail (FT), which are coloured in red, green, blue and magenta, respectively. Stripes indicate 
different outburst stages, see Table[T] The inserts show the pulse profiles at different times in the 2-3.7 and 10-24 keV energy bands (blue and red histograms, 
respectively). 



ing this interpretati on). These component s have been identified 
in other AMSP too (iFalanga et al. 2005a bli lciierlinski & Poutaneril 
I2OO5I : iFalanga et al-l bOOTf T The spectral shapes are very similar 
for individual objec ts over the course of the outbursts (see e.g. 
Gilfanov et ai] 1 19981) as well as for different objects jPoutanenl 
20061) . The accretion disc signatures have been also revealed with 
the XMM-Newton data either spectroscopic ally by the presence of 
the so fter component as in XTE J175 1-305 dGierlinski & PoutanenI 
12005 ^ or by redu ction of the pulse vari ability amplitude below ~2 
keV as in J 1808 jPatruno et d]|2009ah . Also the spectral features 
observed around 6-7 keV in J 1808 are a clear signature of the iron 
line produced by reflection from the neutral matter, presumably the 
accretion disc. 

Understanding of the physical nature of the spectral com- 
ponents has important implications for the correct interpretation 
of the pulse profiles and particularly their strong energy de- 
pendence. Phase-resolved spectroscopy re veals that the two ma- 
jor components do not va ry in phase dGierlinski et"ai] |2002| : 
iGierl inski & Poutanen 200?), resulting in prominent soft time lags 
(i.e. hard photons arriving earlier), which have a steep energy de- 
pendence up to 7 keV, the energy where the contribution of the 
blackbody becomes negligible. The corresponding pulse profiles 
associated with the components also are significantly different, 
with the harder photons showing more harmonic content. A natural 
interpretation of these phenomena is related to the different angular 
pattern of the blackbody and the shock emission (PG03). 

The variability of the pulse profiles in the course of the out- 
burst also gives us a clue to the origin of the X-ray emission. J 1808 
demonstrates remarkably similar evolution of the profiles during 
its outbursts in 1998, 2002, 2005, and 2008 (H08: lHartman et all 



I2OO9I) . During most of the outburst (at SD stage, see Fig.[T), at high 
flux level, the profile is very stable and is nearly sinusoidal with a 
low harmonic content. The profiles observed during the SD stages 
of the 1998 and 2002 outbursts are almost identical (see profiles 
2 and 4 in fig. 3 in H08). Also in 2005, profiles 4A and 4B that 
appeared during a large fraction of the outburst are rather similar 
to those in the SD stages of the 1998 and 2002 outbursts. Even 
in 2008, the pulsar had a si milar profile at high flux period (see 
profiles 1 and 2 in fig. 1 of iHartman et al.|[2009D . Such profile is 
consistent with being produced by only one hotspot as discussed 
and modelled by PG03. The pulse stability allows to obtain a high 
photon statistics and use the average pulse profile to get constraints 
on the neutron star mass-radius relation and the equation of state 
(see PG03 and lLeahv et"ai]|2008l) . 

However, there are clear deviations from this profile both at 
high and low flux levels. In the peak of the outburst in 2002 and 
2005, the profile has a dip in the middle of the broad maximum 
(see profiles for the P3~P5 period in Fig. [!)■ On the other hand, at 
low fluxes, the pulse either has a clearly double-peaked profile or a 
significant skewness opposite to that observed at SD stage (H08). 
This variability also results in jumps in the phases of the harmonics 
(timing noise) with the most dramatic example being a glitch-like 
feature in the phase of the fundamental in 2002 and 2005 outbursts 
dSurderi et al.l200 6: H08). Such variations seriously complicate the 
study of the spin evolution of the pulsar. In 1998 this kind of jumps 
have not been detected as the RX TE data do not cover w ell enough 
the low flux periods (see fig. 1 in iGilfanov et al.|[l993 or fig. 1 in 
H08). 

There are plenty of reasons why the pulse profile might change 
(see |Poutanenll200^ for a review). At high accretion rate, the ab- 
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sorption in the accretion stream might play a role. The effect should 
be largest close to the phase where the flux is large, when the stream 
impact point (i.e. the hotspot) is visible at presumably lowest in- 
clination to the normal. As the accretion rate changes, we expect, 
for example, some variations of the angular emissivity pattern of 
the hotspots, and changes in the s pot area, their shape, as well as 
their position at the stellar surface jLamb et alj2008l : |Patruno et al.l 
l2009c ). as the gas follows different magnetic field lines. Also, the 
visibility of the antipodal spot, hidden at high accretion rates, can 
vary dramatically as the accretion disc retreats, causing a strong 
pulse shape variability. 

Changes in the accretion disc inner radius are expected on 
physical grounds as the magnetospheric radius increases with the 
dropping accretion rate. In addition, there is an indirect evidence 
for increasing inner radius from the drop of the kilohertz quasi- 

J ieriodic o scillation (QPO) frequency during the 2002 outburst 
van Straa ten et al. 2005), if interpreted as a signature of Keplerian 
rotation. 

In the present work we track the changes in spectral and tim- 
ing characteristics through the 2002 outburst of J1808 with the aim 
to determine the variations in the geometry of the system. We first 
analyse the phase-averaged spectra and their evolution. We study 
the variations of the apparent hotspot area and the amplitude of 
the reflection features (which are signatures of cool material in the 
vicinity of X-ray source) during the outburst. We then study the 
pulse profiles, their dependence on energy, the corresponding time 
lags, and finally the phase-resolved spectra. Furthermore we esti- 
mate the physical size of the hotspot and put constraints on the 
geometry, in particular the inner disc radius and the displacement 
of the magnetic dipole from the rotational axis. Using two alterna- 
tive methods, we determine the magnetic moment of the neutron 
star. And finally, we introduce a simple model of the pulsar with 
two antipodal spots and a variable inner accretion disc radius and 
compare the predictions of the model to the data. 



2 OBSERVATIONS 

We focus on the data obtained during the most data-rich outburst of 
J 1808 happened in 2002. The observations are made by the RXTE 
during MJD 52562-52604 (October 15 - November 26) and be- 
long to the ObsID 70080. For the analysis, we used the standard 
HEADAS 6.1 package and the CALDB. We use the data taken by 
RXTE/FCA (3-20 keV) and HEXTE (25-200 ke V); 1 per cent sys- 
tematics has been added to the PCA spectra (see ljahoda et alj2006l 
for a complete review on PCA calibration). To keep the calibration 
uniform throughout our dataset, we used the data from PCA units 
2 and 3 only. 

The outburst can be divided into four periods. The peak (P) 
stage (MJD 52562.1-52563.5; 2002 October 15-16) marks the 
very beginning of the observations, when the 3-20 keV flux (cor- 
rected for absorption) is relatively constant at~ 2.2 x 10~^ erg 
cm~^ s~^ and fast (up to ~ 700 Hz) QPOs are present. The slow 
decay (SD) stage (MJD 52564.2-52574.3; October 17-27) shows 
QPOs twice as slow and the exponential decrease of the flux. At 
the rapid drop (RD) stage (MJD 52574.8-52579.3; October 27 - 
November 1), the flux falls even faster, and the secondary maxi- 
mum appears in pulse profiles after October 29. During the flaring 
tail (FT), the last stage (MJD 52579.4-52602.8; November 1-24), 
the flux was rising and fading with the period of ~ 5 days till the 
source become undetectable. The power density spectrum at this 
stage is dominated by a strong 1 Hz QPO dvan Straaten et al .120051 ; 



Table 1. Data groupings for spectral and timing analysis. 



MJD interval 


Group code 


Outburst stage 


52562.13 


-52563.55 


P 


Peak 


52562.13 


- 52562.46 


PI 


Peak begins 


52562.46 


- 52562.55 


P2 




52563.11 


-52563.14 


P3 




52563.18 


-52563.21 


P4 




52563.25 


-52563.55 


P5 


Peak ends 


52564.25 


- 52574.29 


SD 


Slow decay 


52564.17 


- 52564.50 


SDl 


Slow decay begins 


52564.50 


- 52564.53 


SD2 




52565.09 


- 52565.42 


SD3 




52565.42 


-52565.53 


SD4 




52566.08 


-52566.11 


SD5 




52566.14 


- 52566.45 


SD6 




52567.07 


- 52568.42 


SD7 




52569.04 


- 52569.48 


SD8 




52570.03 


- 52570.28 


SD9 




52570.95 


-52571.39 


SDIO 




52571.94 


-52573.30 


SDll 




52573.33 


- 52574.29 


SD12 


Slow decay ends 


52574.84 


-52579.36 


RD 


Rapid drop 


52574.84 


- 52577.40 


RDl 




52577.40 


- 52579.36 


RD2 




52579.40 


- 52602.79 


FT 


Flaring tail 


52579.40 


-52583.30 


FTl 


Flaring tail begins 


52583.75 


- 52588.44 


FT2 




52588.67 


- 52592.92 


FT3 




52593.02 


- 52602.79 


FT4 


Flaring tail ends 



IPatruno et al.l2009bh . The outburst lightcurve and the sample pulse 
profiles are shown on Fig.[T](see also figs 2 and 3 in H08); note the 
shape differences between various energies. 

While the photon count is rather high at the beginning of the 
outburst, at its end it is necessary to co-add many individual ob- 
servations to obtain the pulse profiles and spectra with reasonably 
small errors. Therefore, we group the data (based on similar count 
rate and pulse shape) as shown in Table[T] For pulse extraction, we 
used the pulsar ephemeris obtained by H08. 

The spectral analysis (both phase-resolved and phase- 
averaged) wa s done using the XSPEC 11.2 spectral package 
( lAmauj[l996h . All uncertainties correspond to a 90 per cent confi- 
dence interval. 



3 SPECTRAL ANALYSIS 

3.1 Evolution of the spectral shape 

The comparison of the spectral shapes at different observation dates 
may help to determine the contribution of various emitting compo- 
nents. Thoughout the paper, we plot the so-called unfolded spec- 
tra, which show the assumed underlying model multiplied by the 
ratio of the observed data to the spectral model folded with the 
detector's response matrix. On Fig. |2] we present the RXTE/PCA 
spectra (unfolded using powerlaw model with spectral index 2.0) 
for several moments of the outburst. In the brightest observations 
at the peak stage we find the flux excess below 4 keV; this fea- 
ture is much less noticeable in the later data. Most probably the 
soft excess is the accretion disc emission. This interpretation is 
supported by the reduced variability at the pulsar frequency be- 
low 2 keV in the XMM-Newton data obtained during the 2008 out- 
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Figure 2. The unfolded spectra for RXTE/PCA range. Vertical lines indicate 
the boundaries for spectral components, from left to right: extra flux below 4 
keV, 6.4 keV Fe line, blackbody hotspot emission below 10 keV and thermal 
Comptonization above 10 keV. The outburst stages coloured as in Fig.[T]are 
indicated at the right. 
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burst jPatruno et alj2009ch . A similar comp onent was found also in 
the XM M-Newton data on XTE J175 1-305 jGierlinski & Poutaned 
l2005h . As the disc temperature is below 0.5 keV in these cases, it 
is not surprising why RXTE sees the excess only at highest fluxes, 
when the accretion disc is expected to be hottest. 

In the 3-10 keV range, the spectrum does not resemble a 
powerlaw and is likely the blackbody-like emission from the neu- 
tron star surface (the similar results were obtained for other AM - 
SPs, see 'Gierliiiski, Done & Barret"2002'; 'palanga et al]|2005dlbl 
iGierlinski & Poutanen 2005; Palanga et al. 2007). Above 10 keV, 
we observe a powerlaw-like emission with the spectral index 
slightly changing in time. The hard X-ray data from RXTEfHEXTE 
(see Section [T2] l displays a cutoff at ~ 100 keV, which is consistent 
with thermal Comptonization in ~ 50 keV electron gas. A bump in 
the spectrum around 6-7 keV and the hardening above 10 keV are 
clear signatures of fluorescent iron line at ~6.4 keV and Compton 
reflection of the underlying continuum from the cold material. The 
amplitude of the corresponding residuals reduces as the flux drops. 

3.2 Phase-averaged spectra 

In this section, we describe the results of fitting the phase-averaged 
spectra with various models. All our models include interstellar ab- 
sorption, with the hydrogen column density fixed at the Galactic 
absorption value of 1.13x10^^ cm~^ (obtained from the HE ADAS 
tool nil using the pulsar's coor dinates). The distance to the s ource 
is assumed to be D = 3.5 kpc jGallowav & Cumming|l200^ . Be- 
cause the relative normalization of the PCA and HEXTE instru- 
ments is uncertain, we allowed this to be an addition free parameter 
in all spectral fits. 

3.2.1 Powerlaw-based models 

A powerlaw model with an exponential cutoff (XSPEC model 
CUTOFFPL), described by the photon index F and cutoff en- 
ergy i?cut, produces rather bad fits. Adding a gaussian iron 
line at 6.4 keV and a blackbody component (BBODYRAD), 
described by the temperature Tbb and normalization K — 



Figure 3. The best-fitting parameters for the blackbody and the cutoff pow- 
erlaw plus reflection model. Evolution of (a) the photon index F and black- 
body temperature together with the flux in 3-20 keV band connected 
for absorption (shown by circles) in units of 10^" erg cm"'^ s~^, (b) 
blackbody apparent emitting area, and (c) reflection amplitude. The vertical 
dashed lines separate various outburst stages. 



[(iioo/km)/(D/10 kpc)]^, significantly improves the fits giving 
X^/doi ~ 1. The apparent residuals remain in the low energy 
channels below 4 keV, at the earlier stages of the outburst (groups 
PI - SD4). We model this feature with another blackbody (repre- 
senting the accretion disc) with the temperature of 0.2 keV and free 
normalization. The resulting fits show that the temperature of the 
main blackbody component decreases from about 0.9 to 0.7 keV 
and the apparent area E — itR%^ also decreases from about 30 
to 10 km^ as the outburst progresses. At the same time the under- 
lying powerlaw softens from F ~ 1.5 to 1.9 and the cutoff en- 
ergy increases from about 40 to 80 keV, and during the FT stage 
it becomes unconstrained. The amplitude of the additional (disc) 
blackbody decreases with time, but it is badly constrained as this 
component is mostly outside of the RXTE energy band. 

The residuals relative the powerlaw fit seen in Fig. |2] im- 
ply the presence of the Compton reflection and an iron line in 
the spectrum. Consequently, we add Co mpton reflection to the 
cutof f powerlaw model (PEXRAV model, iMagdziarz & Zdziarskil 
1 19951) . with the additional fitting parameters being the reflection 
amplitude 5R — ^/(2n) (where Q is the solid angle covered 
by the cold reflector as viewed from the isotropic X-ray source) 
and inclination which we fix at i = 60°. The inclination pa- 
rameter affects only weakly the s hape of the Compton reflec- 
tion continuum in the X-ray ban d i Magdziarz & Zdziarskil [l993 : 



IPoutanen. Nagendra & Svenssod [199^ and does not affect any 
other spectral components. Therefore, any other value of i would 
be also acceptable. 

The limited statistics does not allow to simultaneously con- 
strain 5R and the cutoff energy for the spectra taken only during 
short time intervals. We have fitted the averaged spectrum for the 
SD phase and obtained Ecut = 63^iq keV; hense we assume 
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Scut = 65 keV in all spectra. Our results are shown in Fig.[3l We 
see that the underlying powerlaw index is almost constant F ~ 1.8. 
The apparent blackbody area E reaches the maximum of ~30 km^ 
(corresponding to the radius of only 3 km) at the slow decay stage 
and decreases at the later outburst stages. The reflection amplitude 
is clearly larger at the peak and drops rapidly during the SD stage, 
this explains the apparent softening of the powerlaw seen in the 
model without reflection. The absolute amplitude of reflection de- 
pends on assumed inclination, scaling roughly as oc 1/ cos i. 

3.2.2 Thermal Comptonization model 

The physical origin of the powerlaw-like spectrum with a cutoff 
is most probably related to thermal Comptonization. Because this 
component is pulsating, it is natural to assume that it is produced 
in the accretion shock above the neutron star surface. Although the 
detailed shock structure can be rather complicated, one often ap- 
proximates it by a single temperature plane-parallel slab. To sim- 
ulate the Comp tonized emission, we use the COMPPS model of 
IPoutanen & Sven sson ( 1996.) in slab geometry. We set the inclina- 
tion i = 60° following iGierlihski et alj ( |2002|) . As in powerlaw- 
based models of Sect. 13.2.11 inclination is not constrained from 
the spectral data alone and a particular choice of i will affect only 
the resulting reflection fraction, which scales as oc 1/ cosi. The 
fitting parameters are the temperature and the Thomson optical 
depth of the Comptonizing electron slab Tc and r and the tem- 
perature and the effective emitting area of the blackbody seed pho- 
tons for Comptonization Tsecd and Ehard- Reflection is described 
by the amplitude K and the iron 6.4 keV line (DISKLINE) by the 
normalization. The statistical limitations of the data and the en- 
ergy resolution of PCA do not allow us to fit the inner disc ra- 
dius (that controls relativistic smearing of the line and reflection), 
therefore we fix it (for spectral fits only) on lOrg (statistically ac- 
ceptable for all observations, — 2GM/(? is the Schwarzschild 
radius)Qand assume the radial profile of the emissivity cic . The 
heated surface around the shock is modeled by an additional black- 
body emission component BBODYRAD (see also Gierlinski et al,l 
2002; Pou tanen & Gierlinski 2003; Gierlinski & Poutanen 2005; 
Falanga et al.ll2005bl , l2007h of temperature Tbb and area Esoft ■ To 
account for the excess flux below 4 keV in groups PI - SD4, we 
used an additional (disc) blackbody component with the tempera- 
ture of 0.2 keV and free normalization. 

On the physical grounds, we expect that the blackbody emit- 
ting area Esoft is similar to the area covered by the Comptonizing 
slab Ehard. Attempts to fit these normalizations independently lead 
to unreasonably la rge effective emitting areas of one of the com- 
ponents (see also Gierlinski & Poutanenll2005l sect. 4.5). We have 
fitted the spectrum of group SD, which has the best statistic among 
our dataset, with both Esoft and Ehard being free parameters, and 
obtained Esoft = 751^5 km^ and Ehard = SOlfs km^. Fixing 
the ratio Esoft /Ehard = 2, gives a good fit and results in a reason- 
able area of Ehard = 37 ± 7 km^. We then assumed equal areas 
Esoft = Ehard = E, and were also able to get good fits for all 
data sets and obtained reasonable values for the areas E ~ 30-110 
km^, which imply the apparent spot radius of Roc ~3~6 km. These 
examples clearly demonstrate that on the quantitative level the re- 
sults are affected by the assumed relations between the COMPPS 

^ Here we assume fixed inner disc radius, because of the limited spectral 
data. Later in the paper we show that the variations of the reflection ampli- 
tude and the pulse profiles require the disc inner radius to change. 
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Figure 4. Broad-band spectra of J 1808 and the best-fitting thermal 
Comptonization-based model, (a) Unfolded observed spectrum at the peak 
stage. Solid, dotted, dot-dashed, long dashed and short-dashed curves 
represent the total model spectrum, thermal Comptonization continuum 
COMPPS, blackbody component, reflection with the 6.4 keV line, and an 
additional soft blackbody below 4 keV, respectively, (b) Unfolded observed 
spectra at the P, SD and FT stages. Solid curves represent the best-fltting 
thermal Comptonization model from Sect. l3.2.2l 



and BBODYRAD normalizations. In the following we assume equal 
areas and keep in mind that the systematic error for the spot area is 
about a factor of two, which translates to a 50 per cent uncertainty 
in Rao ■ We also note that to obtain the actual spot size at the stellar 
surface, the apparent areas should be corrected for the inclination 
and gravitational light bending effect, as discussed in Sect. 15. II 

The spectral fits to the broad-band spectra averaged over the 
outburst stages and the contribution of individual spectral compo- 
nents are shown on Fig. |4] We see that the overall spectral shape 
does not see m to vary much durin g the whole outburst as was no- 
tices also by iGilfanov et"ai] ( IT998I) for the 1998 outburst. How- 
ever, the fitting parameters do change significantly (see Fig. |5]l- 
We clearly see that the seed photon temperature Tsocd as well as 
the temperature of the additional blackbody are decreasing with 
time, which is natural as the total luminosity drops. The area E 
also seems to decrease. This trend is independent of the assumed 
ratio Esoft /Ehard- Such a behaviour is expected when the magne- 
tospheric radius increases at lower accretion rate. The values of E 
in the Comptonization model are significantly larger than those in 
the powerlaw model fits (cf. Fig. O, because in the Coinptoniza- 
tion model only a fraction of the blackbody seed photons escape 
through the Comptonizing slab and the Comptonized continuum 
has a low-energy cutoff unlike the powerlaw which continues to 
lower energies without a break. Despite the large uncertainties, the 
reflection 5R shows the decreasing amplitude which is consistent 
with the trend obtained with PEXRAV model. The Thomson optical 
depth and the electron temperature of the Comptonizing slab are 
rather constant, r ~ 1.0-1. 2 and Tc ~40-45 keV, which is consis- 
tent with the stability of the underlying powerlaw in the PEXRAV 
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Figure 5. The evolution of the spectral parameters for the Comptonization model of Sect. 13.2.21 and the coirelations between them. The flux corrected for 
absorption is in 3-20 keV band in units of 10^'-* erg cm~^ s~^. The points are various outburst stages are coloured as in Fig.[T] 



model. The decreasing of the reflection together with the drop of 
the QPO frequencies (if related to the Keplerian frequencies) is 
consistent with the increasing inner disc radius in the course of the 
outburst. 



4 TEMPORAL PROPERTIES 

4.1 Pulse profiles and their energy dependence 

The pulse profiles of J 1808 evolve substantially during the 2002 
outburst (H08). During the peak stage, the energy-averaged pulse 
is close to sinusoidal, but skewed to right. The small secondary 
minimum is noticeable at energies above 10 keV. During the SD 
stage, the pulse is very stable and nearly symmetric, with a slightly 
faster rise than decay. In the end of RDl stage (MJD 52576, Octo- 
ber 29), the secondary maximum (signature of the second hotspot) 
becomes pronounced. At stage RD2 (after MJD 52577, October 
30), the pulse becomes skewed to the left, and the secondary max- 
imum shifts to the rising part of the main peak. On MJD 52581 
(November 3), the pulse has a clear double-peak profile leaving no 
doubts that the second hotspot is visible. At the FT stage, the pulse 
continues to be skewed to the left. The strength of the first overtone 
(i.e. at double pulsar frequency) clearly anticorrelates with the flux 
at later stages as was pointed out by H08. 

We are also interested in the energy dependence of the pulse 
profiles as this provides us with the clues of the origin of pulse 
variability. The evolution during the outburst of the pulse profiles 
in two energy intervals 2-3.7 keV and 10-24 keV is presented in 



Fig. [T] These energy bands are of particular interest, because the 
soft one has the largest contribution of the blackbody component, 
while the hard one contains only the Comptonized continuum. For 
most of the observations, the pulse shapes can well be described 
just by two harmonics: 

F{<f>) = F{1 + ai cos[2-K{ct> - ct>i)] + aa cos[47r((?!. - 4>2)]}. (1) 

We have extracted pulse profiles in various energy intervals for ev- 
ery satellite orbit and fitted them with equation lHJl. The best-fitting 
(positive) relative amplitudes ai, a2 and the phases 01, 02 for the 
fundamental and first overtone for the two energy intervals are pre- 
sented in Fig.|6] We also fitted pulse profiles averaged over the out- 
burst stages in four different energy band (see left panels of Fig.|7]l. 
We see that the profiles depend strongly on energy. During the P 
and SD stages, the relative amplitude of the first overtone is larger 
at hard energies (above 10 keV). The rms drops with increasing 
energy only slightly. At the later stages, the rms at energies above 
10 keV is 2-3 times smaller than that below 4 keV and the oscilla- 
tion amplitude at the fundamental is much more pronounced at soft 
energies comparing to the hard ones (see also Fig.|6]l. 

We also notice that the phase of the overtone behaves differ- 
ently below 4 keV and above 10 keV. There is a clear increas- 
ing trend at low energies, while the phase is consistent with be- 
ing constant at higher energies (compare lower panels in Fig. |6]l. 
This serves as a warning against using such phase-coi mecting solu- 
tions to determine the pulsar frequency evolution (cf. lBurderi et al.l 
|2006|) . Clearly the pulse profile energy dependence and variability 
introduces a strong bias. 
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Figure 6. Fourier amplitudes and pfiases of tlie best-fitting curves (T) to the per-orbit pulse profiles in (a) 2-3.7 keV and (b) 10-24 keV energy bands. Top 
panels: amplitudes of the fundamental ai (squares) and the first overtone 02 (circles). Bottom panels: phases of the fundamental (pi (squares) and the overtone 
4>2 (circles). The vertical dashed lines separate various outburst stages. The typical errors are shown for some points. 



4.2 Time lags 

The pulse profile dependence on energy can also be quantified by 
studying the time lags. Although they contain less information than 
the pulse shape, it is customary to study them as the function of 
energy. In J 1808 the la gs are so ft, i.e. the pulse peaks at softer en- 
ergies at a later phase jCui et al. 1998; Gierliriski et al. 20C^. We 
compute the time lags by fitting the pulse profiles at a given energy 
bin with equation Q and finding the phase difference relative to 
the reference energy^ The 90 per cent confidence interval for one 
parameter is estimated from Ax^ = 2.71. For the four outburst 
stages, the lags relative to the 2-3.3 keV band at the fundamen- 
tal frequency and the first overtone computed from the best-fitting 
phases (pi and (j>2 are presented in the right panels of Fig.|7]by cir- 
cles and triangles, respectively. The absolute value of the lags is 
increasing from 3 up to ~10 keV, and saturates there. Such a be- 
haviour is simi lar to what is observed in the J 1808 d ata during the 
1998 outburst jCui et alJI 19981: iGierlihski et al.ll2002h as well as in 
other sources (see e.g. lGallowavetal.'2002';' Gierliiiski & PoutanenI 
|2005|) . except IGR J0029 1-1-5934, where the lags seem to change the 
trend around 15 keV l lGallowav et alJl2005l ; lFalanga et alj|2005bh . 



^ Alternatively, we can use the standard discrete Fourier transform of the 
pulse profile to compute the relative amplitudes and phases: 

rt 

a^^27.ik4>k = 2/Arph ^ Xje^''^^^''' , (2) 

where n is the number of phase bins, Xj (j = 1 , . . . , n) is the number of 
counts in the jth bin, TVpjj = X]j=i the total number of counts, and 
k = 1,2 coiTespond to the fundamental and the first overtone, respectively. 
The amplitudes and phases obtained by the two methods become similar in 
the limit of small errors and only when the profile can be well represented 
by equation (T). 



The lags at the fundamental increase when the ac cretion rate drops 
(see also iHartman. Watts & Chakrabartvll2008bh . They are nearly 
zero at the peak of the outburst and reach —300 fis at the latter out- 
burst stages. The lags at the first overtone seem to follow similar 
trends, being consistent with zero at the peak of the outburst and 
reaching —200 /is latter on. 

An alternative measure of the lag is the shift with energy of the 
maximum of the best-fitting curve (TJ. To estimate the uncertainties 
of the lag value, we perturb the fitting parameters within 90 per 
cent confidence interval and find the maximal displacements of the 
curve maximum. The results are shown by squares in Fig.|7] These 
lags in general follow the energy dependence obtained from the 
fundamental, except at the peak stage, where they are significantly 
larger. This difference originates in the existence of the absorption- 
like feature at zero phase in the pulse profile at high energies and 
the corresponding shift of the maximum to an early phase. 

The lags have their origin most probably in difference of the 
emissivity patterns of blackbody and Comptonization emiss ion, 
affect ed, in addition, by the fast stellar rotation jGierliilski et al.l 
I2OO2I : PG03). This explains not only why the lag change dra- 
matically below 7 keV, where the blackbody's contribution 
varies significantly, but also the shape of the pulse profile as 
well as its energy dependence. On the other hand, models in- 
volving multip l y Compton down- (or up -) scattering (see e.g. 
ICui et aljri998l : iFalanga & Titarchukll2007l) pay attention only to 
the lags ignoring the pulse profiles. A specific time-lag model of 
iFalanga & TitarcliuS ( l2007h suffers also from other problems. The 
soft lags there are explained by multiply scattering of the hotspot's 
hard radiation in the accretion disc. However, the hard photons 
above 10 keV are actually mostly reflected (by single Thomson 
scattering), while the photons at a few keV are immediately ab- 
sorbed (by photoelectric absorption) because the disc is rather cool. 
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Figure 7. Pulse profiles (left panels) in the 2-3.3, 4.5^.9. 6.5-7.7 and 15.2-22.3 keV energy bands (solid, dotted, dashed and dash-dotted lines, respectively) 
as shown by their best-fitting expression (T). Time lags (right panels) as a function of energy relative to the 2-3.3 energy band for various outburst stages. The 
time lags obtained from the pulse maximum, and the lags for the fundamental and the first overtone are shown by black squares, red circles, and blue triangles, 
respectively. 



Thus, the role of multiply scattering is negligible in any case and 
the lags cannot possibly be produced this way. 



Table 2. Harmonic fits by expression (T) to the phase-resolved apparent 



4.3 Phase-resolved spectroscopy 

The phase-resolved spectroscop y of J 1808 for its 1998 outburst has 
be en performed by[G ierliiiski etaD ( |2002|) and for XTE 175 1-305 
by iGierlinski & Poutaneni (i2005h . It has been concluded, that the 
energy dependence of the pulse profiles and the soft time lags can 
be explained with a simple model where only the normalizations of 
the Comptonized tail and the blackbody vary. The pulse shapes of 
these two components are different with the Compton component 
having stronger harmonic content and peaking at an earlier phase. 
Such profiles in their turn can be reproduced in a physical model 
proposed by PG03, where the angular dependence of the emissiv- 
ity of the spectral components is different. In addition to normal- 
ization, one would expect variations of the blackbody temperature 
caused by Doppler shift, but this is a small effect at the level of 1-2 
per cent (PG03). 

We have generated the phas e-resolved spect r a for the P, SD 
and FT outburst stages. Following lGierlmski et al.l ( l2002h . we used 
the corresponding thermal Comptonization model fits of the phase- 
averaged spectra (presented in Section [3.2.2t as a reference. We 
fixed all the model parameters and fitted the normalizations of 
the blackbody and thermal Comptonization (denoted as i^soft and 
J^"hard, respectively). These normalizations can be directly related 
to the apparent emitting areas (measured in km^) as Esoftihard ~ 
vrD^ifsoftihard (whcrc D — 0.35 is the source distance in units of 
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10 kpc). The results are shown on Fig. [8] The parameters of the fits 
with harmonic function ([T} are given in Table |2] 

At the peak stage, the blackbody variations are consistent with 
a sine-wave, while the Comptonized component shows a dip at zero 
phase, which is reflected in the energy-resolved pulse profiles at 
high energies (left panel, Fig.|7]l. This dip (probably associated with 
absorption in the accretion column) causes the shift of the maxi- 
mum to earlier phase. The profiles du ring the SD stage are almost 
identical to those observed in 1998 dGierlmski et al.ll2002h . with 
the hard tail being less variable and having a stronger harmonic. 
It arrives earlier than the blackbody, resulting in the soft lags. At 
the FT stage, the blackbody shows enormous variability by a fac- 
tor of 3, while the rms variability of the Comptonized tail is still 
only ~5 per cent. Although the absolute value for the blackbody 
emitting area and the amplitude of its variation is model-dependent 
(because of our assumption of the equal phase- averaged areas and 
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Figure 8. Results of the phase-resolved spectral analysis. The best-fitting parameters, except the normalizations of the Comptonization and blackbody com- 
ponents, are frozen at the values obtained for the phase-averaged spectra (groups P, SD and FT). Left panels: effective emitting area of the two components 
(Egoft and Sjiard. blue circles and red crosses, respectively) together with the fits by expression (T) and parameters given in Table|2] Right panels; Sgoft 
versus Shard - The outburst stages are indicated on the panels. 



because the ~0.6 keV blackbody contributes little to the RXTE en- 
ergy band above ~3 keV), it is clear that the blackbody is much 
stronger variable than the tail. These profile shapes are reflected 
in strong variations of the energy-resolved pulse profiles shown in 
Fig.|7l with the low-energy pulses (where the blackbody contributes 
more) having larger rms and the high-energy ones looking exactly 
as the Comptonized tail. The nature of such strong variations of the 
blackbody is not clear. 



5 THEORETICAL IMPLICATIONS 
5.1 Spot size 

The actual spot radius at the neutron star surface can be estimated 
from the apparent spot area S = vri?^ correctin g for the light 
bending and geometry. Using iBeloborodovl ( l2002h analytical ex- 
pression for light bending, PG03 obtained the observed bolometric 
flux produced by a circular blackbody spot of angular radius p (vis- 
ible at all phases) as a function of phase </> of a slowly rotating 
pulsar: 



■A, (3) 



F{(l)) = {l-uf^TvRi sinp^Q + wtan^l + f/cos 

where Jo is the intrinsic radiation intensity, R, is the neutron star 
radius, u = rs/R* and 



U 



u 1 sm I sm f 



It + (1 — u) cos i cos ( 



(4) 



and i is the observer inclination and 6 is the magnetic inclination, 
i.e. the colatitude of the spot centre. The phase-averaged value of 



F{(j}), corresponding to the expression Jsj without the last term in 
brackets, is related to the apparent radius as 

n2 



77f 



D2 ' 



where the observed intensity is I = Io{l — u) 
parent size of the spot we get 

/ 2 P 

Roc, ~ R* sin p IQ + u tan — 



(5) 

Thus for the ap- 



(6) 



which can be reduced to the biquadratic equation for sin(p/2) with 
the only physical solution: 

-1/2 



. p Roc 

sm— — — 

2 i^,^/2 



Q+\ Q 



U COS t COS t 



(7) 



This relation breaks down when parts of the spot are eclipsed, and 
for the homogeneously bright star we have R, = Roo/V^ ~ 
(obviously, then no pulsations can be observed). 

If the apparent size is significantly smaller the stellar radius, 
there exists a simple relation between the physical and the apparent 
sizes (PG03): 

-1/2 



pR* — Roc Q 



(8) 



The smallest possible spot radius is close to the apparent radius as 
(5=1 for i = and 6 — 0. As the oscillation amplitude is small, 
we expect that 9 is rather small, thus the reasonable upper limit on 
the spot size is reached at i = 90° (taking 9 = 0): 



pR, = Ro 



(9) 
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Figure 9. Contour plot of constant angular size (in radians) of the emitting 
spot p for the apparent spot radius Rod = 5 km. Solid curves correspond 
to the neutron star of M = 1.4Mq and radius ij* = 12 km. The actual 
radius of the spot reaches minimum of 5.1 km at j = and 9 = and 
maximum of 8.7 km at i = 90° and d = 90°. The dashed curves corre- 
spond to hi = 1.7M0 and = 11 km. In that case, the actual spot size 
is between 5.1 and 7.6 km. 



5.2 Oscillation amplitude and constraints on the geometry 

Some constraints on the system geometry can be obtained from the 
observed oscillation amplitude. For a blackbody emitting spot at a 
slowly rotating star, the relative amplitude a\ of the fundamental is 
(see PG03 and equation l3l): 

" Q + iitan2(p/2)' ^^^^ 

Comparing this expression to the observed oscillation amplitudes, 
we can constrain the geometrical parameters i and 6. The spot an- 
gular size p is estimated from equation ([7} and substituted to equa- 
tion JlOb . The resulting constraints on i and 6 are shown in Fig. 
llOf a) as the contours of constant amplitude a\ of 5 and 10 per cent. 
We see that for most probable inclinations, the displacement of the 
spot centre from the rotational axis is just 3°-7°. These results 
only we akly depend on the assumed stellar compactness and stellar 
rotation jPoutanen & Beloborodovll2006) . A stronger effect comes 
from the deviation of the emissivity pattern from the blackbody, as 
the Comptonization tail is expected to have a broader beaming pat- 
tern a nd therefore smaller amplitude (PG03,.Viironen & Poutanenl 
|2004) . The same observed ai requir es a slightly larger 9 for a given 
i (see lGierhhski & Poutanei]|2005l) , but still 61 < 10°. Inclination 
smaller than about 40° (and thus possibly larger 6) can be ruled out 
on the bas ic of the X-ray ( 'Chakrabarty & Morgan 1998) and op- 
tical tGiles. Hill & Greenhill 1999 : Homer et al. .2001 . : Wang et al.l 
l200lUDelove et alj2008h orbital modulations. 



Using equation we compute the dependence of p on i and 6 
for a given R^o- The results for R^o ~ 5 km (corresponding to the 
apparent area of 80 km^) and two stellar compactnesses are shown 
in Fig. |9] The actual spot radius pR, is about 5-8 km depending 
on the inclination and compactness. The effect of rapid rotation 
changes the results only slightly. However, the apparent spot size 
obtained from the fitting of the spectra of J 1808 is model depen- 
dent. The powerlaw-based model (Fig. [3} gives an apparent radius 
Roo between 3 and 1 km (for the assumed distance of D = 3.5 
kpc). The radius estimated from the thermal Comptonization mod- 
els (Fig.|5} is significantly larger, varying from about 5.6 to 3.6 km 
(corresponding to E =100-40 km^) during the outburst. The main 
reason for this discrepancy is that the powerlaw continues to lower 
energies without a cutoff, while thermal Comptonization turns over 
at energies corresponding to seed photons. In addition, for Thom- 
son optical depth of the order unity in the Comptonizing slab, only 
about 1/3 of the seed photons can escape directly, the rest being 
scattered and Comptonized to higher energies. Thus on the physi- 
cal grounds we prefer the larger apparent size. 

The apparent area is still uncertain within a factor of two as 
we discuss in Sect. 13.2.21 because of the unknown ratio of the areas 
of the blackbody and Comptonized components. The correspond- 
ing 50 per cent uncertainty in Roc means that the actual spot sizes 
of course are also uncertain by about the same amount, as it is clear 
from equation {Sj. An additional complication comes from the de- 
viation of the spectrum from the blackbody. The actual spot size can 
be larger by the colour correction factor /coi = Tcoi/Tcff , which 
describes the shift of the spectral peak relative to the blackbody. 
However, for the atmospheres heated from above, this correction 
should not play a significant role (see discussion in PG03). In that 
case, the spot size during the whole outburst is smaller than the 
stellar radius for the typical neutron star masses and radii. 



5.3 Visibility of the antipodal spot and the inner disc radius 

The simplicity of the sine-like pulse profiles observed during the 
SD stage can be used to argue that the antipodal spot is not visible 
at this stage. The appearance of the antipodal spot later in the out- 
burst is natural as the accretion disc is expected to recede from the 
neutron star opening the view of the lower stellar hemisphere. If we 
associate the appearance of the secondary peak in the pulse profile 
(on 2002 October 29) with the appearance of the antipodal spot, we 
can estimate the disc radius at this specific moment. This then can 
be used to get independent constraints on the magnetic field. The 
results depend somewhat on the spot size. We first assume Roo ~ 5 
km (corresponding to apparent area of 80 km^, see Fig.|5]l, and then 
vary it by 50 per cent within the uncertainty range. 

Let us consider the parameter space where the antipodal spot 
can be eclipsed by the disc. As parameters we take the inclination i 
and the spot centroid colatitude 6. We take the spot size p given by 
equation Q and shown in Fig. |9] and numerically compute photon 
trajectories from the spot elements towards the observer at every 
pulsar phase and check whether it crosses the disc (see Appendix 
let . This allows us to estimate the maximally allowed inner disc 
radius Rin, which fully blocks the spot. The contours of Rm at the 
plane i-9 for Roo ~ 5 km are shown in Fig. [TOl a). We see that 
at large inclinations i ~ 80°, the disc fully covers the antipodal 
spot even when the inner radius is rather large Rin ~ 40 km. At 
i ~ 50°, the disc has to extend to < 20 km to cover the spot. The 
disc extending to the corotation radius 




blocks the antipodal spot from any observer at i > 67° (for AI = 
1.4M0 and/?, = 12km). 

If we assume a smaller spot Roo ~ 3 km, the inner disc radius 
needed to fully block the spot increases by about 20 per cent com- 
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Figure 10. Constraints on inclination and colatitude of the emitting spot, (a) The solid (nearly vertical) curves show the contours of constant oscillation 
amplitude ai (given by equation IIQI ) of 5 and 10 per cent expected from a blackbody spot at a neutron star of mass M = 1.4Mq and radius Rt = 12 
km. The spot size corresponding to Roo = 5 km is computed using equation 0, see Fig. [9] The dashed curves correspond to the upper limit on the inner 
disc radius (in km) that is required to fully block the view of the antipodal emitting s pot. The corresponding results for the neutron star of M = 1.7Mq and 
Rf = 11 km are shown by dotted curves. The upper limit on the inclination of 82° teildsten & Chakrabartvll200lh is shown by the horizontal line, (b) The 
iimer disc radii required to block the antipodal spot (M = 1.4Mq and ij» = 12 km) for various apparent spot sizes Rao = 3 (solid curves), 5 (dashed), and 
7 km (dotted). 



pared to the iioo = 5 km case (compare solid to dashed curves in 
Fig.llOb). On the other hand, increasing the apparent spot size to 7 
km, necessarily decreases the disc radius by 20 per cent (see dotted 
contours). The stellar compactness influences the results very little, 
for example, taking a higher mass of M = 1.7Mq and i?* = 11 
km, lead to the decrease in the disc radius by only ~ 1 km (compare 
dashed and solid curves in Fig.llOb). 



If the accretion is centrifugally inhibited at Rin > Rcoi, 
the clear signatures of the antipodal spot during the FT stage, 
give us an upper limit on the inclination i < 64° , 67° , 73° for 
Roo —3,5,7 km, respectively (for M = 1.4Mq), which is lower 
than the constraint i < 82° obtained f rom the absence of binary 
eclipses (Bildsten & Chakrabartvll200lh . For a more massive star 
AI — 1.7Mq, the corresponding limits on the inclinations in- 
crease by 3 degrees. If, however, the inclination actually is > 75°, 
this w ould imply that accret i on proceeds even when Rj^ > J?cor 
(see ISBr uit& Taam 19931; iRappaport. Fregeau & SpruitI |2004 



iKluzniak & Rappaport.2007r) 



The constraints on the oscillation amplitudes (see Sect. 15.21 
and solid curves in Fig. 1 10b ) reduce dramatically the allowed re- 
gion in the (i, 6) plane. Combining this together with the fact that 
the antipodal spot has just appeared to the view, gives directly the 
estimation on the inner disc radius at that moment. The results are 
presented in Table[3]for various inclinations and two sets of stellar 
compactness. We see that the largest uncertainty comes from the 
unknown inclination. The compactness affects the results at less 
than 10 per cent level, while a factor of 2 error in the apparent spot 
area (i.e. 50 per cent error in Rao) results in only a 20 per cent error 
in Rin- 



Table 3. Estimation of the inner disc radius and stellar magnetic moment ii 
from the apparent spot area at 2002 October 29. 
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5.4 Constrains on stellar magnetic field 

Having the information about the inner disc radius allows us to es- 
timate the stellar magnetic moment. We need to know the accretion 
rate, which can be obtained from the observed bolometric flux F: 



L = 47rLiV/ang = vMc^ 



(12) 



where L is the total luminosity, rj — 1 — a/1 — u is the accretion 
efficiency in Schwarzschild metric, and /ang is the factor correcting 
for the angular anisotropy of radiation and general relativity effects. 
It depends on the inclination, stellar compactness and spin, as well 
as position of the emitting spot and its emission pattern (see Ap- 
pendix |A] for details of calculations). We approximate the angular 
distribution of radiation from the spot by a law 



/, (a) — Io{l + h cos a) , 



(13) 
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where a is the angle relative to the stellar normal and h is the 
anisotropy parameter The value of h = Q corresponds to the 
blackbody-like emission pattern and the negative h correspond to 
the Comptonized emission from an optically thin slab (see PG03; 
Iviironen & PoutanenI |2004) . For inclinations in the range 50°- 
70° and a circular spot, /ang varies in a rather narrow range 1.3±0.2 
(see Fig. lAlt . It depends very little on the assumed spot size, when 
Rao varies from 3 to 7 km, this factor changes only by 2-3 per cent. 

The inner disc radius depends on the accretion rate M and 
scales with the classical Alfven radius as 



7?in = fcA(2GM)-'''^M-2''V'''^ 



(14) 



where /i is the neutron star magnetic moment. The coefficient fcA is 
rather uncertain, with the numerical simulations giving /ca = 1/2 
jLong. Romanova & Lovelac3l2005l) . Thus, we finally can express 
the magnetic dipole moment as 



0.56 X k 
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(15) 



where ^25 = m/10^^ G cm^. 

The antipodal spot appears after 2002 October 29 at the flux 
level in the 3-20 keV band of F3-20 ~ 0.4 x 
s~^ and the corresponding bolometric flux F ~ 0.8 x 10" 



lO"** erg cm"^ 



erg 

1 



cm~ s~ and luminosity of about L — 1.1 x 10 erg s 
Substituting this to equation dlSt . we get the magnetic dipole mo- 
ments, which are presented in Table [3] for two sets of stellar pa- 
rameters and various inclinations. Rather conservative limits are 
/^25 ~ (9 ± 5) (for anisotropy factor /ang = 1.3 ± 0.2), 

with the largest uncertainty coming from the unknown inclination 
and parameter ^a. The estimated surface magnetic dipole field is 
about Bo ~ (0.8±0.5) x 10* k'jj^'^ G (assuming neutron star radii 
10-15 km), which is in excellent agreement with the value obtained 
recently from the pulsar long-term spin-down rate in the quiescence 
(H08), if fcA '-^ 1. A similar value (close to the upper end) can be 
obtained (see lOilfanov et al1ll998l) assuming that a sharp break (at 
2002 October 28) in the pulsar lig ht curve is associated with the 
onset of the propeller effect (Illarion ov & Sunvaevlll97^ . 

On the other hand, an independent knowledge of the stellar 
magnetic field provides constraints on parameter ^a. Taking ^25 = 
5 ± 3 as obtained by HOB, together with our estimations of the 
inner disc radius from Table [3] we get very conservative limits (for 
inchnations 50° < i < 70°) 



0.4 < fcA < 2.5. 



(16) 



These limits include all systematic and statistical uncertainties, in- 
cluding uncertainties in inclination, stellar mass, anisotropy fac- 
tor, spot size and magnetic moment. As all theoretical models pre- 
dict fcA < 2^^^ ~ 1.1 (see discu ssion in iPsaltis & Chakrabartvl 
I1999I : [iCluzniak & RapDaDortll2007h , this limits the inner disc ra- 
dius to J?in ^31 km and gives some constraints on the inclination 
i < 75°. 



5.5 Varying reflection 

Variation of the reflection amplitude 5R during the outburst and its 
correlation with the flux (see Figs[3]and[5) are clear signatures of 
the changes in the geometry happening when the accretion rate de- 
creases. At the peak stage, our thermal Comptonization fits require 
reflection to be 5R ~ 0.6-0.8 and powerlaw-based models give 5R 




Flux.,_5o (10"" erg cm"'^ s"') 



Figure 11. Reflection amplitude 5ft = f2/(27r) as a function of the flux in 
3-20 keV band. The data points shown by diamonds and stars (with rep- 
resentative error bars) are for the fits using therinal Comptonization model 
of Section [3.2.2l and powerlaw-based model of Section [3.2. II respectively. 
The set of three curves give the theoretical dependence expected for three 
different geometries of the emitting spot: (from top to the bottom) a narrow 
ring with the outer angular radius pmax defined by the dipole formula (TT) 
and the inner angular radius Pmin = O.Qpmax; a ring with inner angular ra- 
dius Pmin = Pmax/2; and a circle of radius pmax centred at the rotational 
pole. The solid curves are for the black body emitting spot (anisotropy pa- 
rameter h = 0), and the dashed curves correspond to the angular pattern 
of the Comptonized emission with h = —0.7, with the magnetic mo- 
ments of /i25 = 6 and 7, respectively, and A'a = 1- Neutron star mass 
is M = 1.4M0 and inclination i = 60°. 



0.4-0.6. The reflection substantially decreases during the SD phase 
and is not constrained at the FT stage. Although the actual reflec- 
tion amplitude is model dependent (depends on the model of the 
underlying continuum and inclination), the trend is not. The dra- 
matic changes in 5R can only be explained if the geometry of the 
accretion flow varies significantly and the inner radius of the ac- 
cretion disc in the peak of the outburst is close to the neutron star 
surface. Assuming i?in in the range 12-15 km at the peak flux level 
of -f3_20 ~ 2 X 10~^ erg cm~^ s~^, we can get a rough estimate of 
the stellar magnetic moment using equation l |15l l. Taking anisotropy 
factor /ang = 1.3 ± 0.2 (see Fig. lAlt and stellar masses in 1.4- 
1.7Mq range, we get ^25 ~ (6±2)fc^^'''*, where the largest uncer- 
tainty comes from a strong dependence on i?in. Such ^ is consistent 
with the estimation from the antipodal spot eclipse in Table [3] 

For illustration, we simulate the expected dependence of the 
reflection amplitude on flux. For various h, we compute the flux 
F as a function of M at a given inclination i (as described in Ap- 
pendix |A]l- The reflection amplitude 5R of course depends on the 
inner disc radius 7?in (see details of calculations in Appendix |B](, 
which varies with M according to equation l ll4b . Both flux and 
reflection depend on the geometry of the emitting spot. As the con- 
straints obtained in Section [5!2l implv a rather small displacement of 
the spot centroid from the rotational axis, for estimating the reflec- 
tion, we can neglect it and assume that magnetic dipole is aligned 
with the rotational axis. We consider three geometrical arrange- 
ments: a full circle centred at the rotational pole, a ring-shaped spot 
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Figure 12. (a) Pulse profiles from two antipodal spots with the inner disc radius varying from 24 to 30 km (from bottom to top). At iJin = 24 km, the 
view of the antipodal spot is fully blocked by the disc, while for = 30 km the whole spot is visible, (b) The amplitudes and (c) the Fourier phases of 
the fundamental and first overtone of the pulse profiles, computed using equation 0, as a function of Ri^- Parameters of the simulations: neutron star mass 
M = 1.4Mq and radius /?* = 2.5rs = 10.3 km, inclination i = 65°, circular blackbody spots of angular radius p =30° centred at a colatitude 6 = 16°. 



with the inner angular radius equal half of the outer radius and a 
very narrow ring. In all cases, we take the maximum spot extend 
from the rotation axis /Omax given by the dipole formula: 



Sm Pinax = 



(17) 



The theoretical reflection-flux dependence are shown in Fig. 
[TT] together with the data. We see that the circular spot does not 
reproduce a sharp dependence of 5R on flux. On the other hand, a 
narrow ring is more consistent with the data, as in that case at high 
fluxes, when the inner disc is close the star, the emission region 
is closer to the equator giving a strong increase in the reflection 
fraction. A sharp increase of reflection and its high value at high 
flux es, can also be e xplained if the inner disc is puffed up (see 
e.g. iLong et al]|2005h because of the interaction with the magne- 
tosphere. 



5.6 Pulse profile variations witii receding disc 

As we argued above, the decreasing area of the emitting spot E and 
of the reflection amplitude 5R when the flux is dropping are strong 
evidences of the receding disc. This is natural as the magnetosphere 
radius is expected to increase. Let us now look at the pulse profile 
variation in the course of the outburst and compare to the expected 
behaviour. 

When the accretion rate is high, the disc should be close to 
the star and may fully blocks the view of the antipodal spot as we 
discussed in Section [531 The accretion flow itself can be optically 
thick and emission from the hotspot can be partially absorbed by 
the flow. Using the mass conservation equation M — iirR^fVp 
(where / ~ 0.1 is a fraction of the stellar surface covered by the 
accretion stream, which corresponds to the spot area of about 100 
km^) and taking accretion velocity ofV^ c/3, we get the Thom- 
son optical depth through the accretion stream at a typical distance 
Rt from the star: 



TT 



pR, 



(Jt 



0.3 



M 



0.1 



1016gS-l / 



(18) 



where ctt is the Thomson cross-section, rrip is the proton mass. At 
the peak of the outburst, the bolometric flux is ~ 4 x 10~^ erg 
cm~^ s~^, and luminosity ~ 6 x 10^^ erg s~^, which gives the 
accretion rate of about 2 x 10^® g s~^ and tt ~ 0.6. Thus, we 



see that the peak of the outburst the flow is marginally optically 
thick and can produce dips in the light curves blocking the hotspot, 
while at lower fluxes corresponding to most of the outburst, the 
flow is transparent. It is possible that this effect is the cause of the 
observed non-sinusoidal pulse shape at the peak of the outburst (see 
the inset for P3-P5 in Fig. [T] and upper left panel in Fig. [TJ. The 
effect is clearly largest at hard energies, which might imply that the 
hard Comptonized radiation is coming from a smaller area (which 
is easier to block) that the softer blackbody radiation. 

As we discussed in the Introduction, the pulse profile during 
SD stage is very stable and is similar to the profiles observed during 
other outbursts at similar flux level. The profiles are well described 
by a one-spot model, implying that the disc is still rather close to 
the star blocking the view of the antipodal spot. The dependence of 
the profiles on energy is consistent with the being produced by the 
different emissivity pattern of the blackbody and Comptonized ra- 
diation modified by Doppler effect (PG03). The pulse stability im- 
plies that the position of the spot at the neutron star surface does not 
change significantly. During the SD stage, the flux drops by a factor 
of 3, which corresponds to an increase of the inner disc radius by 40 
per cent and corresponding increase of the outer spot angular radius 
by 20 per cent according to equations J14t and ( |17| (. However, these 
variations seem to have little impact on the pulse profile, which is 
expected if the centroid of the spot does not move and the spot size 
is smaller than the stellar radius (Poutanen & Beloborodov 2006). 

With the dropping accretion rate, the disc moves sufficiently 
far so that the observer can see at least part of the antipodal pole. 
We thus expect that the pulse profiles change correspondently. Sig- 
nificant variations in the pulse profile are observed from the middle 
of the RD phase. These changes are also reflected in sharp jumps 
in the phase of the fundamental (see Fig.|6l. This motivates us to 
model theoretically variations of the pulse profiles with the varying 
inner disc radius. 

As an illustration, we consider two circular antipodal spots 
emitting as a blackbody and displaced from the rotational axis by 
angle 6. The light curve produced by each s pot is computed follow- 
ing techniques described in PG03 and Poutanen & Beloborodov! 
(2006), accounting for Doppler boosting, time delays, and gravi- 
tational light bending. Now we also account for the effect of ab- 
sorption by the disc of varying inner radius Rin- We compute the 
trajectory for photons emitted by each spot element at every phase 
and check whether it crosses the disc plane at radius smaller or 
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larger than _Rin (see Appendix |C] for details). Fig. ll2r a) shows the 
pulse profile variations as a function of the inner disc radius. No- 
ticeable signatures of the second spot appear when R\n > 27 km 
(for i = 65°), the dramatic change in the pulse profile happens 
when the inner radius changes from 28 to 30 km, as most of the an- 
tipodal spot then appears to our view. These fast changes are also 
reflected in a rather large phase shift A<^i ~ 0.2 of the fundamen- 
tal, but small variations of the phase of the overtone, as shown in 
Fig. [T2l c). This behaviour is qualitatively similar to that observed 
in J1808 (see Figs[6]and|7l(. Thus the phase jumps can be explained 
by changing visibility of the antipodal spot as the accretion disc 
recedes from the neutron star. The timing noise ob served in other 
pulsars (Papitto et al j200H2008l:lRiggio et alj2 008) may also have 
its or i gin in this effect. In addition, spot wandering I Lamb et al., 
l2008l : IPatruno et alj|20093) and variations in the spot size, shape, 
and emissivity pattern might be involved. 

We also show the variation of the amplitudes in Fig. I12f b). 
We see that the fundamental amplitude decreases smoothly as iiin 
grow 13, while the amplitude of the first overtone stays nearly con- 
stant. Such a beha viour is consistent wi th that observed in Fig.[6tb) 
(see also fig. 1 in 'Hartm an et ZlbOOgl) . but is not consistent with 
that seen in Fig. |6{a). We note, however, that in real situation the 
spots are probably not circular and emission is not a blackbody. 



6 SUMMARY 

In this paper, we present the detailed analysis of the 2002 outburst 
data of SAX1808. Below we summarize our findings. 

(i) The phase-averaged spectra of the pulsar are roughly de- 
scribed by a powerlaw with the photon index ~2. However, ad- 
ditional features are significantly detected: iron line at around 6.4 
keV and Compton reflection, as well as a soft blackbody-like com- 
ponent below 5 keV. The dominating powerlaw-like component 
shows a clear cutoff above 50 keV, and the spectral shape is con- 
sistent with being produced by thermal Comptonization in plasma 
of electron temperature ~40 keV. At the peak of the outburst, we 
also detect a soft excess below 3^ keV, which we associate with 
the presence of the accretion disc. 

(ii) The amplitude of Compton reflection is correlated with the 
flux. In spite of the fact, that the actual value for this amplitude is 
model dependent, the trend is not. This behaviour is generally con- 
sistent with the scenario where the disc recedes as the accretion rate 
drops and the magnetospheric radius grows. The observed sharp in- 
crease of the reflection at the highest fluxes is inconsistent with the 
simple circular model of the hard X-ray emitting region, but more 
consistent with a narrow ring. The high reflection fraction and its 
sharp dependence on flux might also imply that the inner disc is 
puffed up because of the interaction with the magnetosphere. In 
any case, the observations require the inner disc radius to be close 
to the stellar surface at the peak of the outburst. 

(iii) The spectral analysis of the phase-averaged spectra also in- 
dicates a strong correlation of the apparent emitting area with the 
flux during the outburst. This result is largely model independent. 



We note that for a slowly rotating star, two antipodal blackbody spots 
which are visible all the time produce no pulsation s whatsoever (see e.g. 
lBeloborodovll2002l:|Poutanen & Beloborodovll2006l) . Thus it is not surpris- 
ing that even for a star rotating at 400 Hz, the amplitude of the fundamental 
decreases when both spots become visible. 



The thermal Comptonization model with equal areas of the black- 
body and Comptonization components give the area dropping from 
about 100 to 40 km^ in the course of the outburst. The actual value 
for the apparent area of the emitting spot depends on the assumed 
ratio of the two components, and has a factor of two uncertainty. 

(iv) We study the variations of the pulse profiles during the out- 
burst and their dependence of photon energy. The profile is very 
stable during the slow decay stage and is similar to those observed 
during other outbursts at similar flux level. However, we observe 
significant variations of the profiles at both high and low fluxes. 

At the peak stage, the deviations are stronger at higher ener- 
gies. We estimate the optical depth through the accretion stream 
and show that at highest fluxes, it is marginally optically thick and 
can absorb part of the radiation at certain pulsar phase resulting in 
the observed dips close to the pulse maximum. The energy depen- 
dence of the effect seems to indicate that the Comptonizing region 
occupies a smaller area compared to the blackbody emitting spot. 

When the flux drops below a certain level, after 2002 October 29, 
a secondary maximum appears in the pulse profiles. This is accom- 
panied by the jump in the phase of the fundamental. We interpret 
these facts in terms of appearance of a secondary antipodal spot. 
We further develop a theoretical model of the pulse profile with the 
varying inner disc radius. We show that varying obscuration of the 
antipodal spot can produce phase jump and pulse profile variability 
similar to those observed. 

(v) Our energy-resolved pulse profile analysis show that the 
pulse shapes at various energies are considerably different. One of 
the quantitative measure of this effect are the time lags, which show 
a complicated variability during the outburst, mostly as a result of 
changing the pulse profile. Relative to the 2-3 keV range, the lags 
are negative (i.e. high-energy emission is coming earlier than the 
softer emission), strongly increasing up to 7-10 keV and saturating 
at higher energies (see Fig. [7}. 

The phase-resolved spectra are well described by two major 
components (blackbody and Comptonized tail) with varying nor- 
malizations, representing the apparent areas of the emission com- 
ponents. In all outburst phases, the blackbody component are lag- 
ging the Comptonized one resulting in soft time lags and explain- 
ing their energy dependence. These facts could be explained by 
the model where the blackbody and Comptonization emissivity pat- 
terns are different. However, the physical nature of the very strong 
blackbody variability observed during the FT stage remains un- 
known. 

(vi) The flux at the moment of significant change in the pulse 
profile, if associated with the appearance of the antipodal spot, can 
be used to put a constraint on the inner disc radius at that moment 
(see TableO- It strongly depend on the unknown inclination i, giv- 
ing radii from about 20 to 35 km, when i varies between 50° and 
70°. The remaining uncertainty, because of the unknown stellar 
compactness and the size of the emitting spot, is about 20 per cent. 

(vii) The information about the flux and the inner disc radius at 
the moment of appearance of the antipodal spot allows us to esti- 
mate the stellar magnetic dipole moment fi ~ (9 ± 5) x lO'^^kjJ^^ 
G cm'^, where fcA is the ratio of the inner disc radius to the Alfven 
radius. The uncertainties here are dominated by the unknown incli- 
nation. A similar value fi ^ (6±2) x 10^^ fc^^^'' G cm^ is obtained 
from the fact that at the peak of the outburst, the inner disc radius is 
in close vicinity of the neutron star as implied by the observed high 
reflection fraction. The corresponding surface magnetic field (for 
NS radii 10-15 km) is Bq ^ (0.8±0.5)10*fcX^''* G. These results 
are in excellent agreement with the values obtained recently from 
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the pulsar spin-down rate (H08) assuming kA ~ 1. To be consis- 
tent with the results of H08, one requires /ca to deviate from unity 
by not more than factor of 2, which give interesting constraints on 
the physics of accretion onto a magnetized star. 

The analysis of the 2002 outburst of J 1808 led us to the follow- 
ing coherent picture of the outburst. At highest observed luminosi- 
ties, the accretion disc is very close to the star. It fully blocks the 
view of the antipodal spot and results in a high reflection fraction. 
The accretion stream at this stage is marginally optically thick and 
absorbs part of the radiation at a phase close to the pulse maximum. 
The behaviour of the reflection seems to indicate a rather narrow 
ring-like emission region and/or puffed up inner disc. As the flux 
drops, during the slow decay stage, the decreasing reflection indi- 
cates that the disc is retreating. The pulse profile is however very 
stable, indicating that the accretion stream is now optically thin and 
the antipodal spot is still blocked from the view by the disc. Dur- 
ing the rapid decay phase, the disc is far enough and the antipodal 
spot opens to the observer's view. This is accompanied by the fast 
changes in the pulse profile and sharp jump in the phase of the fun- 
damental. The outlined scenario is consistent with all the data on 
J 1808 obtained during its five outbursts. 

There are still a number of uncertainties and unanswered ques- 
tions. In our analysis and interpretation of the data, we have made 
somewhat contradicting assumptions. For spectral modelling, we 
have assumed a fixed inner disc radius of lOrg, while the pulse 
variability and changes in the reflection indicate a varying radius. 
This is not a problem, as the obtained reflection fraction is not af- 
fected much by the assumption on the disc radius. For estimations 
of the inner disc radius and stellar magnetic moment from the flux 
at the moment of appearance of the antipodal spot, we assumed a 
circular emitting spot. On the other hand, the reflection-flux depen- 
dence indicates that the correct hotspot geometry might be closer to 
a ring. This difference does not affect the conclusions as long as the 
ring thickness is not much smaller than hal f of its radius, which is 
consistent with the numerica l simulations jRomanova et aljr2004l : 
iKulkami & Romanovall20o3) . Our theoretical calculations of the 
pulse profile variations as a function of the inner radius would not 
be affected much in that case. 

There are other effects not yet accounted in modelling. For ex- 
ample, we have assumed equal areas of the blackbody and Comp- 
tonization components. A dip in the high-energy profile at the peak 
of the outburst might indicate, however, that the high-energy com- 
ponent is produced in a smaller area. This discrepancy does not 
affect our conclusions regarding the decreasing spot area during 
the outburst and the estimations of the inner disc radius as we ac- 
counted for a factor of two uncertainty. However, the calculations 
of the pulse profile will be affected, because the retreating disc 
would affect the pulse profiles at low and high energies in a dif- 
ferent way. This is an interesting problem for further studies. In 
order to reduce the uncertainties still present in the model and to 
get better constraints on geometrical parameters (e.g. inner disc ra- 
dius), it would be useful to know the exact spot shape as a function 
of the magnetic inclination and the inner disc radius. A progress in 
this direction is possible with the detailed MHD s imulations of the 
accretion disc-magnetosphere interaction (see e.g.'Roman ova et aEI 
[2004, ;.Kulkarni & Romanova 2005). A further progress would also 
require a detailed model of the shocked region, where the high- 
energy radiation is produced, together with the radiative transfer 
modelling in order to predict the angular emission pattern from first 
principles. 
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APPENDIX A: FLUX FROM A RING-SHAPED SPOT 

For approximate calculation of the flux and reflection, we neglect 
magnetic inclination, and consider the emitting region in the form 
of a ring symmetric around the stellar rotational axis. The bound- 
aries of the ring are described by colatitudes pmin and pmax. We 
assume that the light from the antipodal ring below the equator is 
blocked by the accretion disc. 

Let us first compute the observed flux from the ring. Con- 
sider a surface element AS = i?JdcosSd(ji at a slowly rotat- 
ing star. We neglect here the effects related to rapid rotation be- 
cause the flux averaged over the ring is affected by them very 
little. The element position is described by the unit vector n = 
(sin 6 cos (j>, sin 9 sin (j), cos 9) that points to it from the star center. 
Let the unit vector along the line of sight be fe = (sini, 0, cosi) 
with i being the inclination. The angle between n and the line of 
sight is denoted by i/; so that 



2.0 



cos tp — k ■ n = cos i cos 9 + sin i sin 9 cos 0. 



(Al) 



The observed flux from this surface element in approxima- 
tion of b eloborodov (2002) for light bending is (see PG03; 
IPoutanen & Beloborodovi,2006) : 
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Figure Al. The anisotropy correction factor as a function of inclination 
computed from equations iA5\ and iA7\ . The solid curves correspond to 
stellar parameters M = 1.4Mq and _R» = 12 km, while the dashed 
curves to M = 1.7Mq and _R» = 11 km. More inclined curves are for 
blackbody-like radiation pattern with h = and almost horizontal curves 
are for h = —0.7. Here a circular spot centred at the rotational pole with 
Roo = 5 km is assumed. 



where the angular distribution of intensity at the stellar surface 
I, (a) is assumed to follow the dependence l ll3t . a is the angle 
of the photon emission relative to the stellar normal related to other 
parameters as 



cos a ~ u + {1 — u) cos — Q + U cos (j), 

with Q and U given by equation 

Integrating l lA2b over the ring surface, we get 

E 



F{t) 



(A3) 



(A4) 



where E = i?,27r(/i2 — Mi) is the ring area, fii = cos pma 
fi2 = cos Pmin, and 

1 



C = u + 
+ h 



'-{ill + ^2) cosi 



+ u{l — u){fii + P2) cosi + i(l — it)^ sin'^ i 



+ /ii(l - ■u)^(3cos^ i - 1) (/ii + P1M2 + (A5) 
6 

which is strictly valid when the whole ring is visible. 

The total emitted luminosity at the stellar surface (for two 
symmetric about equator, ring-shaped emitting regions) is 

L. = 2E 27r ^ /. (q) cos a d cos q = 2E Jo 27^^ + ^^, (A6) 

and the total luminosity at infinity L = (1 — u)L,. We thus get 
the anisotropy correction factor relating the observed flux to the 
luminosity in equation J12b : 



/ang — 



1 I h 

2 3 



(A7) 



(A2) 



47rD2f(i) (I-m)C' 
As an example, let us consider a circular spot (i.e. pmin ~ 0) 
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centred at the rotational pole. For a given apparent spot size Roo , 
we compute the physical spot angular size pmax ~ p from equation 
Q. The dependence of /ang on the inclination is shown in Fig. lAll 
Variations of Roa within 3-7 km, lead to only 2-3 per cent changes 

in /ang- 



APPENDIX B: REFLECTION FROM THE ACCRETION 
DISC 

Let us now estimate the reflection amplitude expected for the emit- 
ting ring. By definition it is given by the ratio of the luminosity 
reflected from the disc to that directly observed: 



r{ip) = 



- cos tp)^ 



1 1/2 



4(1 + COS?/')2 



+ 



sin^ ip 



rs(l — cos tp) 
2(1 + cost/;) ' 



(B9) 



Substituting here ■(/) = i/>d, we obtain the disc radius, where photon 
trajectory crosses the disc plane. 

Calculation of the reflected fraction now involves simple inte- 
gration over the ring area (upper hemisphere) and photon emission 
angles: 



d cos 6 / J, (a) cos a d cos a / d</3 _H"(r[i/)d]— -Rin), (BIO) 



where H is the Heaviside step function and the limits on a and 
are given by conditions l lB4t and l |B5t . 



— ./"ang^^ri 



■A, 



(Bl) 



where Prcfl is the fraction of emitted photons that hit the disc sur- 
face at radii r > Rin- 

Because of the axial symmetry, the surface element in a ring 
is described only by colatitude with the corresponding unit vector 
pointing towards it, n = (sin ^, 0, cos S). The photon direction 
(in the frame related to the element, with 2:-axis along its normal) 
is described by two angles, polar angle a and azimuth ip measured 
from the southern direction of the meridian. At large distances from 
the star the photon direction is given by vector k. For a given a, the 
angle its makes to n can be obtained from equation l |A3l ). Using 
spherical trigonometry it is easy to show that 

(cos ?/j sin 9 + sin cos 9 cos if \ 
sin ip sin ip 1 . (B2) 

cos Tp cos 9 ~ sin ip sin 9 cos if J 

The necessary condition that the trajectory crosses the disc is 



cos Ip cos 9 < sin tp sin 9 cos ip, 



(B3) 



that translates to the limits on a and p. In the stellar polar region 
sin 6* < K. = u/{l — u) (existing if it < 1/2), the limits are 



< (/p < 27r, 
cos p > cos po , 



if < cos a < cos a_ , 

if cosQ- < cos a < cos Qa 



(B4) 



where cos a± = u±(l — it) sin 9 and cos ipo = cot ip cot 9. Thus 
photons emitted close to the surface will always cross the disc plane 
irrespectively of the azimuth, and those emitted close to the zenith 
cos a > cos a+ do not cross it at all. In the intermediate case, 
photons have to be emitted in the southward direction to cross the 
disc plane. Outside the polar region sin 9 > k, the limits are 



cos(^ > cos^o, if < cos a < cosa+. 



(B5) 



Now let us compute the radius where the photon trajectory 
crosses the disc plane. The unit vector along the intersection line of 
the disc and photon trajectory planes is 



d 



{cosip, cos6'siniy5, 0) 
y/l — sin^ 9 sin^ ip 



and thus the angle it makes with k is given by 

sin tp cos 9 + cos tp sin 9 cos ip 



cos — d ■ k 



^/l — sin^ 9 sin^ p 



(B6) 



(B7) 



For the given emission angle a and the impact parameter 



b ^ R,sma/VT^, (B8) 
an approximate photon trajectory is ( lBeloborodovll2002h : 



APPENDIX C: ANTIPODAL SPOT ECLIPSE BY THE 
ACCRETION DISC 

The pulse profiles produced by the spots on a rapidly rotating 
star are computed fo llowing techniques described in PG03 and 
IPoutanen & Beloboro dov ( 2006). However, the effect of the accre- 
tion disc needs to be accounted for, because the photons emitted by 
the lower antipodal spot can be absorbed on the way to the observer 
by the disc. We divide the spot into a number of elements and com- 
pute the light curve from each of it independently. The position 
of each element can be described, as in Appendix |A] by the unit 
vector n = (sin S cos 0, sin ^ sin 0, cos 6*), where now 9 > n/2. 
The angle tp this vector makes with the direction to the observer 
k = (sin i, 0, cos i) is given by equation JAU . 

The element in principle (without the disc) is visible to the 
observer at all phases if cos{9 + i) > — k, not visible at all if 
cos{9 — i) < — K, and visible only when coatp > —Q/U if 
cos(9 + i) < —K < cos{9 — i) (see IPoutanen & Beloborodovl 
l2006l for details). If the element with the given {9, cp) can be vis- 
ible, we substitute tp given by equation JAll ) to the light bending 
formula ([A3} getting the emission angle a and the impact parame- 
ter b using equation dBSI l. which fully define the photon trajectory. 
Its plane intersects with the disc plane along the line described by 
the unit vector 



d = 



(— cos f 



COS t sm y cos < 



' COS i sin (p, 0) 



cos^ i + cos^ 9 — 2 COS i cos 9 cos ip 

The angle it makes with k is 

, , , COS i cos Ip ~ cos 9 

cos ip^ = d ■ k — 



\J cos^ i + cos^ 9 — 2 COS i cos 9 cos ip) 



(CI) 



(C2) 



We find the radius, where trajectory crosses the disc r('i/)d), using 
equation ( |B9t . If it is larger than J?in, the element is invisible. 



